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Abstract 

In this paper we consider the classical and quantum control of 
squeezed states of harmonic oscillators. This provides a method for 
reducing noise below the quantum limit and provides an example of 
the control of under-actuated systems in the stochastic and quantum 
context. We consider also the interaction of a squeezed quantum os- 
cillator with an external heat bath. 



1 Introduction 

In this paper we consider the problem of squeezing of harmonic oscillators 
from the point of view of control theory. Squeezing has been suggested as a 
method for reducing noise in quantum systems below the standard quantum 
limit. This can be achieved by using laser pulses and in that sense may 
be viewed as a quantum control problem, although the classical squeezing 
problem is also of interest. In the latter case one is interested in reducing 
noise induced by random perturbations. 

The quantum control problem has been of great interest recently, see 
for example Brockett and Khaneja [1999], Lloyd [1996] and Warren et. al. 
[1993] and references therein. 

Here we consider squeezing as a control problem in both the classical 
and quantum setting. In the classical case we consider a system subject 
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to thermal noise while in the quantum case we consider a system at zero 
temperature and in the presence of noise. In both cases the control is given 
by an external electromagnetic field and enters the control equations multi- 
plicatively. In this sense the setting is similar to the NMR control problems 
analyzed by Brockett and Khaneja. 

A key feature of squeezing is that it results in a redistribution of uncer- 
tainty between observables. 

In this paper we consider a model for phonon squeezing in solids following 
the work of Garret at. al. [1997], but one can equally well consider the case 
of photons in quantum optics. The control is via a single pulse on a large 
ensemble of oscillators and this sense we are considering under-actuated 
control systems in both the classical and quantum case. 

We also model the effect of dissipation on the classical system and the 
effect of coupling to a heat bath in the quantum setting. This causes the 
squeezing effect to gradually moderate. 



In the classical setting we consider bilinear control systems of the form 



where a^- is a constant matrix (i.e. the free dynamics is linear), bij is linear 
in x, i.e. the control u 3 enters bilinearly, w k is white noise and the state 
space is W l . 

In the quantum context we want to consider a similar equation but 
defined on an appropriate Hilbert space: 



where H s is the Schroedinger operator, the Hj are (linear) input operators, 
u> are functions of time, and the ip is a vector in the Hilbert space. 

3 Classical Squeezing of the Harmonic Oscillator 

In this section we consider classical squeezing of a set of identical coupled 
harmonic oscillators. Denote the position of each oscillator by u l . 



2 The Control Setting 
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The Hamiltonian for the system is of the form: 

i i 

where the oscillators are assumed to have unit mass and pi = u l . 

In order to analyze the system we decompose it into its normal modes. 
Denoting the normal mode coordinates by Q l we thus obtain a system of 
uncoupled harmonic oscillator equations of the form Q l + QfQ 1 = 0. 

The main control mechanism we consider here is squeezing by pulses. 
In this case each oscillator is forced by a pulse at time t = which is 
proportional to its displacement, i.e. we have equations of the form: 

Q l + n?Q* = 2XQ i S(t) (3.2) 

where S(t) is the Dirac delta function and A is a constant which is propor- 
tional to the frequency Q. 
Thus we obtain 

Q\0 + ) = Q l ((T) + 2AQ*(0) 
Thus, if one considers the system subject to white noise, 

Qi + nfqi = 2XQ i 8(t) + oii* , (3.4) 

one sees that while one starts with a spherical equilibrium distribution 
which is invariant in time, after the pulse one has an elliptical distribution 
which rotates in time at twice the harmonic frequency (by the Z2 symmetry 
of the ellipse). (A precise analysis is given below in the course of our treat- 
ment of the quantum mechanical case.) Noise reduction is then achieved by 
viewing the system "stroboscopically" when the noise is low. 

Actually the above is an idealization: in actuality the oscillator should 
be viewed as in equilibrium with a heat bath which dissipates energy. In 
the classical setting one can model this by simple linear dissipation (in the 
quantum setting one has to introduce a heat bath - see below). 

Thus we have a system of the form 

Q* + nfQ i = -TkQ* + Ui(t) + avj (3.5) 

where rji is a dissipation constant and E/j (t) is the control which we can choose 
to be a single pulse or a sequence of pulses. Depending on the dissipation 
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strength an initial squeezing effect will decay away and we need a continual 
sequences of pulse to keep the system in a squeezed state. 

It is worthwhile remarking on the how the control enters in our setting: 
the control is a single pulse applied overall (and in this sense the system is 
under-actuated) while the effect on each (normal mode) of oscillation is to 
apply a pulse proportional to displacement (minus the mean displacement 
which is of course zero for each oscillator). This is effected by the type of 
interaction of the oscillators with the field that the pulse induces. We note 
also that in the full nonlinear setting the mean displacement may not be 
zero and must be taken into account. 

We shall return to the classical squeezing of oscillator by pulses, and in 
particular a computation of mean square displacement, after a discussion of 
the quantum case below. 

We note also that parametric resonance control can achieve similar squeez- 
ing effects in the classical case. In this case we consider oscillator motion in 
the presence of a modulating drive: 

Q-Lo 2 (l + ecos2ojt)Q = (3.6) 

where e parameterizes the strength of the drive. We omit details of this 
approach here. 



4 Squeezing of the Quantum Harmonic Oscillator 

We now turn to the quantum setting. 
Consider the following Hamiltonian 

H = ^ + ^fQ 2 + \5{t)Q\ (4.1) 

which reflects an impulsive change in the spring constant and where oo = 
y K/m, K being the original spring constant. 

The variables P and Q, which are operators in the quantum case, obey 
canonical commutation rules [P, Q] = ih. We can rewrite the above Hamil- 
tonian in terms of creation operators a and defined through 

p =V¥<< ,t -<'>> < 42 » 

with [a, at] = 1. Written in terms of the new variables, the Hamiltonian is 
H = hw(a) a + 1/2) + \5{t){a + a)) 2 . (4.3) 



The ground state of the system, for t 7^ 0, |0), corresponds to the vacuum 
of a, (a|0) = 0), and the excited states are of the form (o^) 2 |0). 

We now want to study the behavior of the system at t > 0, given that the 
system is in its ground state at t < 0. The wave function at t = + is of the 
form \ip(t = + )) = exp(— i\Q 2 )\0), and for longer times the system evolves 
with the "unperturbed" Hamiltonian: \tp(t > 0)) = exp(— iHot)e~ lX Q |0). 
Our first quantity of interest is (ip(t)\Q 2 \ip(t)) = (Q 2 (t)). Let us compute it 
using the general method of coherent states. We find 

(Q 2 (t)) = {Q\e aQ \ae- iut + a)e' iut ) 2 e- iXQ2 \Q), (4.4) 

where we have used the fact that e tHot ae~ tHot = ae~ tujt , which states that 
a) and a respectively destroy and create eigenstates of Hq, and where Q is 
defined in units of ^Jhj \2muj). 

Now we introduce a basis of coherent states \z), which satisfy a\z) = 
z\z), { z \ a ^ = ( z \ z * i an d form an overcomplete set of states: 

1 = ^- [ dz dz*e- zz *\z)(z\. (4.5) 



Inserting (4.5) in (4.4) we find 



(Q 2 (t)) = -L f [dz dz*e~ zz * 

27TZ J J 

(z 2 e~ 2iujt + z* 2 e 2iujt + 2zz* - 1)| (0|e 4A ^ \z) \' 

In order to evaluate the last term we need the position representation of the 
ground state (note that at this point Q is a real number) 

(0\Q) = \e-® 2 / 2 (4.6) 

7T4 

and that of the coherent state 

(Q\ z ) = J_ e -Q 2 /2W2,Q^ 2 /2. (4 . 7) 
7T4 

A simple integration gives 

{0\e iXQ2 \z) = [ dx{0\Q)(Q\z)e iXQ2 (4.8) 



Changing to the variables z = u + iv we have 



— zz 



*Me lXQ2 \z)\ 2 



1 



e -[v 2 + (2X 2 +l)u 2 +2Xuv]/(l+\ 2 ) 



(4.10) 
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VT+X 1 



and 



^ POO /»CO 



tt cos ujt + v sin wt + m> sin 2wt — - 



) 



xe -[« 2 + (2A 2 +l)« 2 +2AH/(l+A 2 ) 

= 1 + 4A 2 sin 2 ut + 2A sin 2ujt 



(4.11) 
(4.12) 



It is interesting to compare this with an ensemble of classical oscillators 
with initial conditions taken from a heat bath. For simplicity let us take 
uj = m = Ub = T = 1 [k,B is Boltzman's constant). An arbitrary oscillator 
will evolve as Q(t) = ucost + vsint, with u and v its initial position and 
velocity. If a pulse is applied at t = of the form treated above: Q(t) = 
u cos t + (v + 2\u) sin t. Now let us average over initial conditions taken from 
a measure given by (a thermal bath) 



It is interesting to note that the two expressions for, respectively, the 
quantum oscillator at zero temperature and the classical oscillator at finite 
temperature, are exactly the same. The general time dependence of the 
variance for a squeezed harmonic oscillator with frequency oj can thus be 
written in the following form 



with eo = hw/2 for the quantum case and eo = ksT for the classical oscillator 
at a temperature T. 

The method of coherent states presented above has the advantage of 
being suitable for calculating other quantities. For example, if the oscil- 
lators are atoms within a solid, the scattering amplitude for an X-ray is 
decreased by a factor (called the Debye- Waller factor - see Ziman [1972]) 
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~ (expikQ(t)), with k the wave- vector of the X-ray. We now ask ourselves 
what is the time evolution of the Debye- Waller factor for a squeezed phonon. 
This means that we need to compute the following expression 

J(A,t) = (0|e^ 2 e( ae " <ut+ate " <wt )e- a « 2 |0) 

= —= , / du dv 

v/e VTTX 1 ^ J 

2u cos ut+2v sin cot- -^ 2 + (2A 2 +i)^+2A^] 
g l+4A 2 sin 2 u;t+2A sin 2u)t 

For the Debye- Waller factor, we obtain the following time dependence 

( e <*Q(*)) = e -k 2 {Q 2 (t)) 

Measurement of the Debye- Waller factor may provide a practical method 
of detecting the squeezing phenomenon experimentally. 
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5 Squeezing and dissipation 

In this section we consider the squeezing of a quantum oscillator coupled 
to a an infinite number of oscillators representing a "heat" bath. We show 
that this causes a decay in the squeezing oscillation for small time and true 
damping in the limit of a continuum of oscillators. This damping effect of 
the heat bath is similar to that analyzed classically in Lamb [1900], Komech 
[1995], Sofer and Weinstein [1999] and Hagerty, Bloch and Weinstein [1999]. 
We stress that we are considering a zero temperature case, and the damping 
effects appear due to a) the coupling of a single variable with a continuum 
of variables and b) an "asymmetry" in the initial conditions. The applied 
pulse on the oscillator generates outgoing waves on the continuum system 
which in turn gives rise to a positive damping (for a detailed discussion of 
negative versus positive damping see Keller and Bonilla [1986]). 

We start with a general formulation, and at the end of this section discuss 
a specific continuum example. 

The Hamiltonian of the system consists of three parts: #0 describing 
the original oscillator: 

pl m*4 2 



the Hamiltonian H e of the environment: 



Pi_ , frmj 2 
2m + 2 9a 



and a linear coupling between the two 

#int = ^ £,aq a qo ■ 



(5.2) 



(5.3) 



Formally, the total Hamiltonian H = H$ + H e + H int can be written in terms 
of its normal mode coordinates X v and P v : 



r p2 



2m 2 " 



(5.4) 



and we will consider a situation in which the initial (before the pulse) wave 
function corresponds to all the modes in the ground state: 



*»=nQ) 1/v "^ /2 " 



(5.5) 



At t = a pulse is applied to the (original) oscillator, the wave function 
immediately after the pulse given by: 



^ (t = + ) = e iXq o^ 

= e iX T, tlv Uo t iUovX l _ l X l ,^^ 



(5.6) 
(5.7) 



where U^ v is the matrix transforming from the original (uncoupled) modes 
to the coupled system (qo = Uq u X u ). 

As in previous sections, we are interested in the fluctuations of the vari- 
ance of qo, given in this case by 



(5.8) 



/IV 



and that we will compute by solving the equation of motion obeyed by the 
correlations (X^X u )(t). Since X^ and X u correspond to harmonic coordi- 
nates, using the quantum mechanical commutation relations we compute 



the equations of motion to be: 



j t (X,X u ) = ±-{{P,X u + P u X,)) 



= ~™ ^l(X,P u )+u; 2 u (X u P,)) 
^{P^Pu) = -{ul + uD^PX + 2mWy„(X fl X u ) . 

Note that the above equations are identical to those of classical harmonic 
oscillators for the quantities Xn(t)X u (t) etc., with initial conditions given 
by the values of the correlations evaluated for the quantum wave function: 



(P„P„)(0+) = ^^ 

+2h A (1 + 5^) — ^ q 

muj^ muj u 

({X^ + P v X li ))(0+) = A\W^U^{— + — ) 

with ql = (q 2 (0-)} = £ a WlJlmuj a . 

Collecting the above equations we obtain 



<?o(*)> = ?o |l + 4A^(t) + -^C(t)S(t)j , (5.9) 
with 

S(t) = V — ^ sin wJ C(t) = V — 9£ cos w„t. 
^— ' raw„ ^— ' mw„ 

All the information of the evolution of the variance is contained in the 
function J(u), the physical interpretation of which is that of a local density 
of states of the oscillator, defined as 

J H = £ -=rr s (» - ( 5 - 10 ) 

— mujii 



from which 



S(t) = J dujj{uj) sin tut, C(t) = J du,J(uj) cos ujt. (5-11) 

Note that J(u>) is a sum over delta functions, giving rise to a superpo- 
sition of oscillations with the frequencies lo v for both S(t) and C{t). In the 
limit of an infinite system, and when the modes are spatially extended over 
all space J(w) becomes a continuous function. In that case the oscillatory 
behavior acquires a damped component, the detailed time dependence be- 
ing given by the frequency spectrum of J{lo). A lorenzian shape for J '(uj) 
will give an exponentially damped oscillation for both S(t) and C(t). As an 
illustration of this point we consider a model for which J(lo) can be com- 
puted explicitly - see the classical analysis in Lamb [1900] Komech [1995]. 
Consider a one-dimensional string coupled to our oscillator. The string is 
described by a "transverse" displacement u(x,t). The classical equations of 
motion of the system are 

Uu(x,t) = c 2 u xx (x,t) 
Md 2 q {t)/dt 2 = -Vq (t)+T[u x (0+,t)-u x (0-,t)} 

q (t) = u(0,t). (5.12) 

The normal modes consist of even and odd (in x) solutions. The odd solu- 
tions do not involve go an d are of the form u q ^ (x,t) = e tcqt smqx, whereas 
the even solutions are of the form u q ^ e (x, t) = e icqt cos(q\x\ + 6 q ), with 8 q 
a phase shift (to be found). The wave vectors q label the normal modes, 
and play the role of the index /i in the above discussion: u>^ = cq, and 
U 2 q = cos 2 (5 q ) (up to a normalization constant) in the present case. Sub- 
stituting this expression in fl5.12| ) we obtain (ujq = V/M) 

tmij= Mc(^|)_ (513) 
from which U 2 = cos 2 5 q is given by 



„2, ,2 



UnO — 2 2,(2 2^2 = (5-14) 



where we have defined a = 2T/Mc. Note that U q represents the transfor- 
mation matrix that has to be normalized and since the frequencies form a 



continuum we normalize U q {u)q) to its integral over uj q . Omitting the index 
q in ojg, we obtain 

,2 



TTI \ 2a U " mUJ Tl \ 

U (") = — „,9. 3. , 7, .9. ^2=— J ( UJ )- ( 5 - 15 ) 



7T a 2 UJ 2 + (u 2 — UJ 2 ) 

Substituting (|5.15|) in (pll) we obtain 



S(t) = e _i * sin Q t, C(t) = e _i * cos fi *, (5.16) 



with 

fin 



k;o (l + [a/wn] 2 ) ^ cos 6/2, (5.17) 
(l + [a/wn] 2 ) ^ sin 5/2, (5.18) 



where 5 = tan -1 a/wn- 

In the realistic limit a -C wo which corresponds to a "weak" coupling to 
the environment) this expressions take the form: S(t) = (h/(mLVo)exp(—Tt/Mc) sinwni, C(t) 
(h/ (m(Jo)exp(- Tt/Mc) cos wnt. 

Note that in this model, and in the limit of weak coupling, the initial 
variance q 2 . of the reference oscillator is unchanged due to the coupling to 
the environment, and is given by q^ = h/2mujQ. Our final result for this 
section is then 



(ll(t)) 

2Xh\ . „ / 2Xh\ . 
mujQ J \ .,, . i 



1 + e -2(T/Afc)t 
2 

mu>o J 



(5.19) 



which reduces simply to ( 4.13| ) in the uncoupled case of T = 0. 

In summary we have shown in this section that the coupling to the 
environment can be included in general, giving rise to dissipation, and that 
the squeezing effect in the presence of dissipation can be computed explicitly 
for the Lamb model. 

Additional details of the analysis here, extensions to the squeezing of a 
nonlinear oscillator, and a treatment of the quantum measurement issue will 
appear in forthcoming publications. 
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